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摘 要 :为 了 研究 开 式 转子 在 低速 起 飞 工 况 下 的 噪声 特性 ,运用 计算 流体 力学 (CFD) 和 计算 气动 声 
学 (CAA) 相 结合 的 方法 ,对 开 式 转子 的 离散 噪声 做 了 数值 模拟 。 采 用 非 线性 谐 波 法 (NLH) 对 开 式 
转子 的 非 定常 流 场 进 行 计算 , 非 定 常 分 量 和 叶 面 定常 载荷 一 同 构建 类 比 声 源 ,最 后 使 用 FW-H 积分 
法 进行 远 场 嗓 声 辐射 计算 ,获取 远 场 多 点 的 骂 声 频谱 结果 。 搭 建 了 开 式 转子 试验 侣 ,使 用 0° ~90。 
远 场 指向 性 阵列 测试 了 指向 角 区 间 内 的 骂 声 特性 结果 ,并 与 数值 模拟 结果 进行 对 比 。 研 究 结果 表 
明 :干涉 噪声 是 开 式 转子 的 主要 噪声 ;关键 频率 处 的 数值 模拟 结果 与 试验 结果 趋势 基本 一 致 ,误差 
在 3dB 以 内 ;所 建立 的 非 线 性 谐 波 法 频 域内 展开 求解 非 定常 扰 动量 的 气动 噪声 数值 模拟 方法 ,可 以 
二 为 开 式 转子 骂 声 的 快速 准确 预测 提供 一 定 的 帮助 。 
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Numerical simulation and experimental measurement of open-rotor 


WEI Kai,CAO Qi,YAN Qun,XU Jian, XUE Dongwen 


( Aviation Technology Key Laboratory of Aeroacoustic and Vibration, 


Aircraft Strength Research Institute,710065 Xian,China) 


Abstract :In order to the study the characteristics of open-rotor noise in low speed, CFD (computational 
jfluid dynamics )/CAA (computing aerodynamic acoustics) hybrid method of the open-rotor tone noise was 
set up. N-S equation is expanded by NLH method in the frequency domain ,and the unsteady components 
are obtained. These unsteady components and steady load are used to construct Lighthill sound source. Fi- 
nally, FW-H integration method is used to calculate far-field noise radiation. The noise experimental plat- 
form is built and test scheme is designed. The noise directivity in the range of 0°-90° is measured by using 
the far-field directivity array surround scale prototype. The characteristic of the noise of open-rotor is ana- 
lyzed by comparing the numerical simulation results and experimental results. The results show that inter- 
action noise is the major component of open-noise. The simulation results are basically consistent with the 
distribution trend of the experimental results at key frequency ,the error is within 3 dB. The numerical sim- 
ulation method of aerodynamic noise based on nonlinear harmonic method is developed to solve unsteady 


disturbance in frequency domain ,which can provide some help for rapid and accurate prediction of open 
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开 式 转子 发 动机 兼 具 涡 桨 发 动机 耗 油 率 低 和 涡 


声 贡献 更 大 。 然 而 随 着 对 噪声 预计 精度 的 进一步 提 


扇 发 动机 适 于 高 速 飞行 的 特点 。 相 比 于 同等 推力 等 
级 的 常规 涡 扇 发 动机 , 开 式 转子 发 动机 耗 油 率 和 
CO, 排放 可 降低 20% ~ 30% ,可 为 最 大 马赫 数 为 
0.70 ~0.85 的 民用 客机 、 军 用 运输 机 提供 动力 1。 
随 着 国际 民航 组 织 对 发 动机 环保 性 要 求 的 进一步 提 
高 ,以 及 更 加 严格 的 适 航 条 例 的 出 台 , 开 式 转子 发 动 
机 将 是 下 一 代 民 用 飞机 潜在 的 动力 备 选 方 案 之 一 。 
开 式 转子 发 动机 有 油耗 低 的 优点 ,但 也 因 其 转子 直 
接 暴露 于 空气 中 ,存在 噪声 大 的 缺点 。 后 者 是 其 走 
向 进一步 成 熟 并 在 商业 化 应 用 进程 中 面临 的 重要 挑 
战 多 一 所 。 开 式 转子 发 动机 对 转 转子 很 大 程度 上 决 
定子 该 型 发 动机 的 飞行 性 能 ,同时 又 是 整个 发 动机 
的 专 要 噪声 源 ; 因此 ,对 商业 运营 和 适 航 条 例 而 言 ， 
开 式 转子 的 噪声 必须 控制 在 一 定 范围 内 。 

CI 目前 ,国内 外 控制 开 式 转子 噪声 的 方法 主要 是 
通 寺 气动 声学 设计 来 加 以 改善 。 该 方向 的 研究 集中 
其 羡 式 转子 气动 噪声 的 产生 及 传播 机 制 、 开 式 转子 
的 噪声 数值 模拟 与 试验 测试 。 文 献 [5] 采 用 计算 流 
体 沪 学 (CFD) 与 计算 气动 声学 (CAA) 相 结合 的 方法 
对 形式 转子 在 起 降 阶段 的 频 域 噪声 进行 了 数值 模 
摄 e 非 定常 流 场 计算 采用 非 定常 雷诺 平均 方法 
(WEANS) ,噪声 计算 采用 基于 Mohring 的 声 类 比方 
法 忆 结 果 表 明 开 式 转子 的 峰值 辐射 噪声 主要 集中 在 
低 镀 范 围 。 文 献 [6] 探 索 了 开 式 转子 的 几何 建 模 方 
法 ,并 基于 URANS 方法 ,针对 巡航 工 况 下 不 同 迎 角 
时 开 式 转子 的 流动 特性 进行 了 数值 模拟 ,未 对 噪声 
进行 分 析 。 结 果 表 明 , 两 级 转子 间 存 在 明显 的 气流 
干涉 现象 , 迎 角 和 两 级 转子 间距 是 影响 开 式 转子 气 
动 特性 的 关键 参数 。 文 献 [7] 基 于 CFD-CAA 的 数 
值 模拟 方法 ,利用 旋转 叶片 的 周期 性 ,采用 单 通道 模 
拟 了 开 式 转子 的 非 定常 流 场 ,提高 了 流 场 计算 的 效 
率 , 采 用 Ffowcs Williams-Hawkings( FW-H) [8] 预测 噪 
声 ,计算 结果 与 试验 结果 吻合 很 好 ,并 且 得 出 结论 ， 
开 式 转子 的 噪声 幅 值 主要 产生 在 转子 叶片 的 通过 频 
率 及 其 谐 波 频率 处 。 文 献 [ 9] 梳理 了 3 个 典型 的 开 
式 转子 模型 噪声 风 洞 试验 ,详细 分 析 了 离散 噪声 和 
宽频 噪声 对 开 式 转子 噪声 的 贡献 程度 。 结 果 表 明 ， 
在 低频 区 域 ,离散 噪声 贡献 大 ; 在 高 频 区 域 ,宽频 噪 


高 ,需要 刻画 更 细节 的 非 定常 流动 ,传统 的 非 定常 模 
拟 技 术 URANS 的 计算 时 间 和 硬件 成 本 昂贵 而 难以 
被 工程 界 所 接受 。 因 此 ,更 高 效 的 非 定 常 建 模 技 术 
一 直 是 CFD 计算 的 研究 热点 。 

非 线 性 谐 波 (nonlinear harmonic , 简称 NLH ) 法 
是 一 种 建 模 思想 相对 简单 ,精度 较 高 ,近年 来 比较 流 
行 的 叶轮 机 械 非 定 常 计算 方法 "1。 该 方法 的 基本 
思想 是 将 流动 分 量 分 解 为 时 均 流 场 和 脉动 流 场 , 然 
后 在 频 域内 求解 确定 性 应 力 项 ,只 需要 对 一 个 叶片 
通道 生成 网 格 ,计算 量 大 幅 降低 ”。 王 雷 等 采 
用 非 线性 谐 波 法 对 对 转 压气 机 的 两 级 转子 的 干涉 现 
象 进 行 了 分 析 。 研 究 表明 非 线性 谐 波 法 可 在 时 均 框 
架 内 有 效 计算 非 定 常 现象 ,3 阶 谐 波 已 能 对 在 计算 成 
本 和 计算 精度 间 获 得 最 佳 折 囊 。 药 晓 江 等 "采用 非 
线性 谐 波 法 和 传统 的 双 时 间 步 法 对 风扇 和 涡轮 分 别 
进行 非 定常 计算 分 析 。 分 析 表 明 在 相同 的 计算 机 硬 
件 条 件 下 ,两 种 方法 达到 相同 的 计算 精度 , 非 线 性 谐 
波 法 可 节省 40% 的 计算 时 间 。 开 式 转子 存在 强烈 的 
级 间 干 涉 特 征 , 目 前 采用 NLH 法 计算 开 式 转子 的 非 
定常 扰动 量 进而 预测 噪声 的 研究 还 比较 少 '™] 。 

本 研究 将 以 某 典 型 开 式 转子 为 研究 对 象 ,基于 
CFD-CAA 的 数值 计算 方法 对 开 式 转子 的 离散 噪声 
进行 数值 模拟 ,采用 非 线 性 谐 波 法 在 频 域 内 展开 求 
解 非 定常 扰动 项 ,与 叶 面 定常 载荷 一 同 构 建 类 比 声 
源 ,最 后 使 用 FW-H 积分 法 进行 远 场 噪 声 辐射 计算 。 
对 比 数值 模拟 与 试验 结果 ,验证 采用 非 线 性 谐 波 法 
在 频 域内 展开 求解 非 定常 扰动 项 的 气动 噪声 数值 模 
拟 方法 的 有 效 性 ,为 开 式 转子 的 噪声 快速 评估 提供 
一 定 的 帮助 。 


1 计算 模型 与 数值 方法 


1.1 计算 模型 


开 式 转子 的 设计 参数 参考 NASA 报告 ,转子 
直径 刀 =0.6 m, 两 级 转子 均 为 8 个 叶片 ,级 间距 离 
0.35 m, 旋 转轴 长 2.4 m,75% 叶 高 处 的 安装 角 均 为 
35°。 开 式 转子 几何 模型 如 图 1 所 示 。 
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图 1 开 式 转子 几何 模型 
Fig.1 Open rotor geometric layout 

流 场 计算 时 ,两 级 转子 各 取 一 个 叶片 组 成 单 通 
道 计 算 域 ,划分 六 面体 结构 化 网 格 ,计算 域 如 图 2 所 
示 。 入 口 边 界 至 桨 中 心 距 离 为 1 了 ,出口 边界 距离 桨 
中 心 3D。 径 向 半径 取 4D, 为 2.4m。 单 通道 网 格 如 
图 3 所 示 。 第 一 层 网 格 高 度 为 10“m, 进 出 口 流 道 
数 均 取 121 ,叶片 区 域 采用 0 型 贴 体 网 格 , 叶 顶 间 阶 
采用 蝶 形 网 格 , 叶 顶 与 叶 根 加 密 网 格 ,采用 3 重 网 格 
技术 ,网 格 总 量 890 万 。 边 界 条 件 总 温 293 ,总 压 
10233 kPa, 
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2 图 2 流 场 计算 域 

C9 Fig.2 The computational domain of flow field 
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图 3 单 通道 流 场 网 格 

Ee Fig.3 Single channel flow field mesh 


加 开 式 转 子 的 噪声 主要 由 厚度 噪声 和 载荷 噪声 组 
成 "1 。 计 算 噪 声 时 ,以 叶片 国 壁面 为 积分 面 ,由 于 流 
场 网 格 尺度 足以 满足 噪声 计算 的 网 格 要 求 , 可 直接 
采用 流 场 计算 时 的 叶片 表面 网 格 ,无 需 局 部 加 密 。 
FW-H 求解 时 直接 重 构 每 个 气动 变量 的 值 ,提取 前 3 
阶 叶片 表面 静 压 作为 输入 ,不 需要 进行 任何 额外 的 
插值 过 程 。 声 学 网 格 如 图 4 所 示 。 

边界 条 件 采 用 周期 性 边界 条 件 ,前 后 排 各 取 一 
个 叶片 组 成 单 通道 。 计 算 域 进口 和 出 口 均 采 用 压力 
边界 条 件 , 总 压 101. 3 kPa , 总 温 293 上, 汕 流 黏度 
0. 1 mm/s 。 轴 向 进 气 ,来 流 无 扰动 。 噪 声 传播 方向 
沿 每 个 网 格 块 指向 外 法 线 方向 。 计 算 远 场 噪声 辐 
射 ,需要 确定 合理 的 噪声 监测 点 ,监测 点 以 声 源 中 
心 -前 级 转子 旋转 中 心 ,为 坐标 原点 ,以 旋转 面 正 前 
方 5 m 处 为 测 点 1,5 m 为 半径 , 画 90" 圆 弧 , 每 间隔 
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10° 为 一 个 测 点 。 


到 4 ”声学 网 格 


Fig.4 Acoustic mesh 
1.2 计算 方法 


本 研究 采用 商业 软件 进行 定常 和 非 定常 数值 模 
拟 。 为 了 适应 开 式 转子 相对 转动 的 流 场 ,定常 数值 
模拟 需要 对 守恒 形式 的 Navier-Stokes 控制 方程 组 '“ 
进行 时 均 化 处 理 , 速 度 相 关 量 均 采 用 相对 值 , 即 可 得 
到 相对 直角 坐标 系 下 的 雷诺 平均 N-S 控制 方程 组 ， 
各 列 向 量 的 具体 表达 式 为 


[Pp | pei Fr 0 
po1 poiw1 +P61 TH 
U=| pw, |, fi;=| poiws +p |, -fy = Ty 
PO3 pOiW3 + pos Ta 
Lpk [| (ph +p) ow; Lg; + wr 
(1) 


式 中 :U 为 守恒 型 变量 ;f 为 离散 的 对 流 项 ;p 为 密 
度 ;E 为 总 能 ;w; 为 速度 项 ;7; 为 应 力 项 。 

定常 与 非 定常 数值 模拟 的 汕 流 模型 均 选 取 S-A 
模型 , 转 静 交接 面 均 采用 一 维 无 反射 混合 面 法 。 

本 研究 采用 非 线 性 谐 波 法 开展 非 定常 数值 模 
拟 ,其 基本 原理 在 于 将 非 定常 流 场 分 解 为 时 均 流 场 
和 扰动 流 场 两 部 分 进行 耦合 求解 。 将 式 (1) 中 的 解 
向 量 世 分解 为 时 间 平 均值 与 春 干 周期 性 扰动 之 和 ， 
如 下 式 所 示 。 

U(r,t) = U(r) + SD U(r,t) (2) 

式 中 :r 为 位 置 矢 量 六 为 时 间 。 每 一 个 周期 性 扰动 
基于 时 间 均 可 分 解 为 N 阶 谐 波 函 数 之 和 , 即 


U(r,t) = 多 Ur)exp(iwt) + 
各 (3) 


U(r)exp(iw) -J 


式 中 ,为 谐 波 阶 数 。 当 =1 时 ,其 频率 等 于 叶片 的 
通过 频率 。 阶 数 越 高 ,对 求解 的 硬件 资源 需求 越 大 。 
为 了 对 时 均 方程 中 的 确定 性 应 力 项 进行 求解 ， 
还 需要 引入 额外 的 非 定常 输 运 方程 。 详 细 推导 见 文 
献 L[11] 。 对 于 每 一 个 周期 性 扰动 量 四 均 可 得 到 一 
个 相应 的 扰动 输 运 方程 ,而 对 于 周期 性 扰动 量 的 每 
一 阶 谐 波 均 可 得 到 其 相应 的 谐 波 方程 为 
0, =- S78 + SI +00 (4) 
式 中 :4 为 网 格 单元 ;了 为 对 流 项 ; 玉 , 为 苛 性 项 ;s 
为 面积 ; 为 源 项 。 
通过 求解 上 述 谐 波 方程 ,确定 性 应 力 可 以 直接 
根据 谐 波 函 数值 计算 得 到 。 对 于 扰动 量 广 和 8& ,确定 


(5) 
EB, |/; | 和 | g, | 为 扰动 量 第 阶 谐 波 幅 值 的 实 
部 站 本 研究 的 开 式 转子 离散 噪声 ,只 需 计算 厚度 品 
堆 和 载荷 噪声 。 提 取 前 3 阶 谐 波 的 压力 和 速度 , 带 
太 吏 进 后 的 FW-H 方程 , 则 

四 


六 | fpov,C7) Te. ds(y) dr + 

_ -oo SCT) (6) 

© | | pr) Sods) dr 

| -oa S(7) 7 

式 惠 :。 为 叶片 的 法 向 速度 ;f 为 叶片 表面 压力 。 
O 通 过 式 (3) 和 式 (5) 可 计算 出 叶片 表面 的 压力 
和 速度 ,再 联 立 式 (6) , 便 可 得 出 开 式 转子 的 远 场 噪 


声 声 压 。 
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2 试验 研究 


2.1 开 式 转子 试验 台 


为 了 在 地 面 声学 环境 中 开展 模拟 开 式 转子 发 动 
机 的 噪声 试验 研究 ,需要 模拟 开 式 转子 在 起 飞 初 始 
阶段 ,发 动机 面临 低 转 速 和 低 来 流速 度 的 特征 。 设 
计 并 人 研制 了 7 开 式 转子 试验 台 , 如 图 5 所 示 。 开 式 转 
子 试验 全 主要 由 转子 系统 、 动 力 驱动 系统 和 控制 系 
统 组 成 。 转 子叶 片 几 何 参 数 如 图 1 所 示 , 与 文献 
[18] 一 致 ,试验 台 设 计 转 速 1 000 r/min。 


ChinaXiv 合 作 期 刊 479 


多 5 半 消 声 室 噪声 测试 照 


Fig.5 Test photo of semi-anechoic chamber 


2.2 试验 方案 


本 试验 在 半 消 声 室 进行 , 开 式 转子 试验 台 位 于 
半 消 声 室 前 部 。 试 验 通过 远 场 阵列 测试 研究 开 式 转 
子 在 0 ~2 000 Hz 范围 内 的 噪声 特性 ,试验 现场 如 
图 5 所 示 ,麦克 风 阵 列 测 点 示意 图 如 图 6 所 示 。 


测 点 11 


测 点 12 \ 测 点 10 测 点 9 
出 


图 6 麦克 风 阵 列 测 点 示意 图 
Fig.6 Test points arrangement 
如 图 6 所 示 , 远 场 麦克 风 阵 列 和 开 式 转子 轴线 
在 相同 水 平面 上 , 以 后 排 转 子 中 心 为 圆心 ,半径 为 
5 m, 在 0° ~110" 的 圆 绝 上 均 布 12 支 自由 场 传 声 器 ， 


传声器 戴 防 风 球 , 即 可 获得 远 场 噪声 数据 。 试 验 时 
同时 连续 采集 12 支 传声器 30 s 的 声 压 数据 ,采样 频 
率 设 为 16 384 Hz ,频率 分 辨 率 为 0.5 Hz。 试 验 分 别 
测量 开 式 转子 在 500 r/min .600 r/min 和 700 r/min 
的 噪声 特性 。 


3 结果 与 分 析 


3.1 数值 模拟 结果 分 析 


数值 模拟 了 开 式 转子 在 2 000 r/min 的 厚度 品 
声 .载荷 噪声 和 总 噪声 ,图 7 和 图 8 分 别 展示 了 测 点 


480 


1 ~5 的 厚度 噪声 、 载 荷 噪声 频谱 图 。 从 图 6 和 图 7 
中 相同 频率 、 相 同 观 测 点 上 可 以 看 到 ,厚度 噪声 比 载 
荷 噪声 小 很 多 。 在 图 7 中 可 以 看 到 , 厚度 噪声 在 
250 Hz 处 为 10 dB , 同 频率 下 的 载荷 噪声 为 15 dB。 
从 图 8 中 可 以 看 出 ,各 离散 频率 处 , 测 点 1 的 声 压 级 
最 大 , 测 点 2、3、4.5 的 声 压 级 依次 减 小 ,说 明 开 式 转 
子 的 远 场 噪声 声 压 级 最 大 的 位 置 在 旋转 平面 正 前 
方 , 即 与 旋转 轴线 呈 0° 位 置 处 。 
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图 7 厚度 噪声 频谱 图 
Fig.7 Thickness noise SPL levels 
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图 8 载荷 噪声 频谱 图 
Fig.8 Loading noise SPL levels 

图 9 展示 了 开 式 转子 在 2 000 r/min 的 总 噪声 ， 
包含 厚度 噪声 和 载荷 噪声 , 与 图 7 和 图 8 对 比分 析 
可 得 载荷 噪声 占 总 噪声 的 绝 大 部 分 分 量 , 厚 度 噪 声 
只 在 低频 区 域 有 少量 的 贡献 。 前 级 转子 的 通过 频率 
记 为 BPF1 ,后 级 转子 的 通过 频率 记 为 BPF2 ,各 转子 
的 谐 波 频率 可 记 为 n x BPF, 其 中 为 谐 波 阶 数 , 取 
正 整 数 。 两 级 转子 各 阶 谐 波 频率 自由 组 合 得 到 的 频 
率 为 开 式 转子 的 干涉 频率 "”。 从 图 9 中 可 以 看 到 声 
压 级 峰值 分 布 在 两 级 转子 的 通过 频率 和 干涉 频率 
处 。 还 可 以 看 出 , 声 压 级 最 大 的 3 个 频率 均 在 干涉 
频率 处 。 对 比 可 知 ,干涉 噪声 声 压 级 较 大 ,是 开 式 转 
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子 的 主要 噪声 。 
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图 9 总 噪声 声 压 级 频谱 
Fig.9 Total noise SPL levels 


3.2 试验 结果 分 析 


试验 开展 了 500 r/min .600 r/min、700 r/min 共 
3 个 转速 工 况 的 远 场 噪声 测试 。 图 10 展示 了 开 式 
转子 700 r/min 时 部 分 测 点 的 声 压 级 频谱 ,在 0 ~ 
1 500 Hz 的 频率 范围 内 , 声 压 级 的 3 次 幅 值 均 出 现在 
通过 频率 及 其 谐 波 干涉 频率 处 。 也 可 以 看 出 ,各 测 
点 的 噪声 频谱 没有 被 转子 通过 频率 及 干涉 频率 所 完 
全 支 配 ,图 中 出 现 了 很 多 窄带 频谱 , 宽频 噪声 明显 。 
这 与 本 次 试验 处 于 低 来 流速 度 、 低 转速 状态 有 关 。 
开 式 转子 在 此 阶段 发 生 了 短暂 的 失速 现象 , 随 着 飞 
行 速度 的 增加 , 既 型 气动 迎 角 和 来 流速 度 会 逐渐 匹 
配 ,升力 也 逐步 增加 。 
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10 ”部 分 测 点 声 压 级 频谱 图 
Fig. 10 SPL levels of some test points 

图 11 展示 了 工 况 700 r/min 时 , 测 点 1 测 得 的 
噪声 频谱 图 ,并 与 同样 条 件 下 的 数值 模拟 结果 进行 
了 对 比 ,主要 对 比 了 通过 频率 ( BPF1 或 BPF2) .干涉 
频率 ( BPF1 + BPF2 )、 及 其 干涉 频率 的 谐 波 (2 x 
BPF1 + BPF2 2 xBPF1 +2 xBPF2)4 个 关键 频率 处 
的 结果 。 可 以 看 出 ,在 关键 频率 处 ,数值 模拟 结果 与 
试验 结果 趋势 基本 一 致 ,误差 不 超过 3 dB ,误差 值 相 
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对 较 小 ,吻合 较 好 。 
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Fig. 11 


SPL comparison of test and simulation at test point 1 
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基于 非 线性 谐 波 法 ,建立 了 开 式 转子 噪声 预测 
模型 ,数值 模拟 了 其 离散 噪声 特性 。 开 展 了 开 式 转 
沪 进 场 噪声 试验 并 与 预测 结果 对 比 , 主要 结论 如 下 。 
COD1) 开 式 转子 的 离散 噪声 主要 由 干涉 噪声 组 成 。 
< 小 噪声 主要 产生 于 两 级 转子 的 通过 频率 及 其 
后 处 。 

加 2 ) 关键 频率 处 ,数值 模拟 的 开 式 转子 离散 噪声 
眶 混 验 值 趋势 基本 一 致 ,误差 在 3 dB 以 内 ,满足 工 
程 精度 要 求 。 

><3 ) 形 成 的 从 时 域 求 解 - 频 域 展开 -构建 声 源 项 - 远 
场 强 射 传播 的 方法 链条 是 有 效 的 ,可 为 后 续 快 速 地 
噪 震 评估 提供 帮助 。 
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